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Abstract 

We consider the ADM formalism of the Brans-Dicke theory and propose an an- 
isotropic extension of the theory by introducing five free parameters. We find that the 
resulting theory reveals many interesting aspects which are not present in the original 
BD theory. We first discuss the ghost instability and strong coupling problems which 
are present in the gravity theory without the full diffeomorphism symmetry and show 
that they can be avoided in a region of the parameter space. We also perform the post- 
Newtonian approximation and show that the constraint of the Brans-Dicke parameter 
cjbd being large to be consistent with the solar system observations could be evaded 
in the extended theory. We also discuss that accelerating Universe can be achieved 
without the need of the potential for the Brans-Dicke scalar. 
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1 Introduction 



There are many alternative theories and extensions of the Einstein's general relativity. 
Motivations to consider them are diverse. They range from a classical straightforward 
extension of general relativity like the scalar-tensor theory to quantum gravity such as 
string theory. A comprehensive review can be found in a recent article [1]. The classical 
modification is based on the observational rationale that the predictions of the general 
relativity has been successful only on the solar system scale and has never been tested on 
the cosmological scale. Recently, it has another input from cosmology, that is, to search 
for the theoretical foundation of the accelerating Universe [2] has become one of the most 
fundamental problems in modern cosmology. There are two main avenues. The first one 
is to assume an unknown source of energy which is repulsive in nature referred to as the 
dark energy and various proposals for its origin have been put forward [3]. The other is to 
consider an alternative theory of the Einstein's theory of general relativity and modify the 
gravity in the IR limit. We are still awaiting more accurate observational data to distinguish 
between them. Therefore modification of general relativity is a very important subject and 
is attracting a great deal of interest with extensive research activities being reported (see [1] 
and references therein). 

Among the many alternatives, scalar-tensor theory [1] is especially interesting. In this 
theory, the scalar field enters in a nontrivial manner, specifically through non-minimal 
coupling term and it is based on solid foundation of general relativity. Typical examples 
would be the Brans-Dicke (BD) theory [5] and the gravity with a dilaton field arising for 
instance in the string theory. The scalar field in this type of theory provides many distinctive 
theoretical aspects, for example like induced gravity [6] and explanation of the behavior of 
the universe in the early inflationary stage as well as the late stage |1] in the cosmological 
context. 

Especially, we attempt to modify the BD theory which is known as the prototype of the 
scalar-tensor theory and is one of the simplest alternative to Einstein's theory of general 
relativity (GR). BD theory was firstly designed to properly incorporate the Mach's principle 
into General relativity by replacing the gravitational constant G with a scalar field 
which can vary with space and time. It is well-known that the observational constraints 
on BD theory is restricted by the astronomical tests in the solar system, i.e., the BD 
parameter cjbd > 50000, being obtained from the observable Post-Newtonian parameter 
7obs-l = (2.1±2.3) X 10-5 [318]. 

Modifications of Brans-Dicke theory have been actively pursued. Most notable is the 
introduction of the potential for the Brans-Dicke scalar. For example, in Ref. [Hllin], BD 
theory with a suitable self-interacting potential for the BD scalar field was proposed to 
account for the accelerated expanding of the Universe. In another example, it can be shown 
that f{R) gravity in the metric formalism is equivalent to BD theory with a potential and 
the parameter wbd = [TT] . 

Another alternative is an anisotropic extension. It is motivated by the Horava-Lifshitz 
(HL) gravity [12] which was proposed as a quantum gravity theory. It is based on anisotropic 
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scaling of space and time as a fundamental symmetry, abandoning the Lorentz symmetry 
at short distance. In the IR limit, the theory reduces to GR when the Lorentz violating 
parameter A becomes 1. There has been considerable interest in the theory and various 
aspects have been investigated [13]. Among others, a BD type of generalization of the 
Horava-Lifshitz gravity was shown to be possible within the detailed balance condition 
in Ref. [H], and the resulting theory reduces in the IR limit to the usual BD theory 
with a negative cosmological constant. From the theoretical point of view, breaking of 
the full diffeomorphism invariance coming from the anisotropy puts strong constraint in 
the resulting theory. The breaking causes the theory to have an additional scalar degree 
of freedom, whose behavior could create serious problems, such as a ghost or classical 
instability or the strong coupling problem. On the experimental side, it should pass the 
solar system test which requires wbd has to be big in the original BD theory. 

In this paper, we propose a more straightforward modification of BD theory starting 
from the ADM formalism and discuss the above issues. A close look at the ADM formalism 
reveals that the original BD theory can be anisotropically extended by introducing four 
more free parameters. Depending on the values of these parameters, the theory shows 
many interesting aspects which are not present in the original BD theory. For example, 
the large value of ubd in the original BD theory to be consistent with the solar system 
observations could be evaded in the modified theory. It can be also shown that the ghost 
instability and strong coupling problems which are present in the gravity theory without 
the full diffeomorphism symmetry can be avoided within some range of the free parameters. 
Another distinctive feature is that the accelerating Universe could be achieved without the 
help of some potential contrary to the original BD theory. 

The paper is organized as follows. In Sec. 2, we consider anisotropic BD (aBD) action 
with five parameters which reduces to the ordinary (isotropic) BD action for particular 
choices of the four parameters. In Sec. 3, we consider a quadratic action through the 
perturbation analysis of the aBD gravity in order to investigate the pathological behaviors 
of scalar graviton mode,. In Sec. 4, we perform the perturbation analysis up to cubic order 
and check whether the strong coupling problem can be cured to this order. In Sec. 5 we 
study the observational constraint on the aBD theory and compare with the experimental 
results. The conclusion and discussions are given in Section 6. 

2 Anisotropic Brans-Dicke action 

In order to construct an anisotropic Brans-Dicke gravity with additional parameters, let 
us first consider an isotropic Brans-Dicke gravity in the ADM formalism, whose metric is 
parametrized by 

dsl^^ = -N^dt^ + gij{dx' + N'dt){dx^ + N^dt), (2.1) 
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where is the lapse function, Ni is the shift function, and gij is the three dimensional 
metric. For the ADM metric, Brans-Dicke action is given by [1] 



+<PR - 2g'^ViVj<P - 0JBD<p-'g''Vi'pVj(l> \ , (2.2) 



where (j) and ubd is the Brans-Dicke scalar and Brans-Dicke parameter respectively, gij and 
R are the three dimen 
and 71 take the forms 



R are the three dimensional metric tensor and the Ricci scalar. The extrinsic curvature Kij 



K^, = ^{9^,-^^N,-V,N,), (2.3) 

7r = ^(0-V,0iVO, (2.4) 

where the dot denotes differentiation with respect to t. The first line of (12.20 is the kinetic 
(K) term and second line is the potential (V) term. 

We first notice that the BD gravity fl2.2p preserves the full diffeomorphism symmetry. 
We can extend it to aBD theory by introducing free parameters which explicitly break the 
diffeomorphism invariance and consider the most general action as follows 0: 



where 



5'aBD — S^BD + 'S'^BD) (2-5) 

^ _ j dt£xN^ {0 {KijK'^ - XK"^) - 2r]iKTi + ^10-^^} , (2.6) 

S^BYi = j dtd^xNy^ {(PR - 2r/2ViV> - a;20"'Vi0V>) . (2.7) 

In this action, the parameters A, ?7i^2; 1^1,2 are dimensionless constants. Note that when 
At^I, ?7i7^1, 7727^1, and Ui ^ uJi, the action ( 12.50 does not have the full diffeomorphism 
invariance, but is invariant under the foliation-preserving diffeomorphism: 

= -C(t,x), 5t = -/(t), 

^Qij = diC^gjk + djC^gik + C^dkgij + fgij, 

6N, = ^,CN, + C'^,N, + Cg^J + fN, + fNi, 

SN = CdjN + fN + fN, 

6<j) = C'5fc0 + /0. (2.8) 



^This extension is based on the fact that each term of 2Ktt, 2g'^^WiWj(j>, (p ^tt"^, and (f> ^ g'^^ W i<f>W j (f> in 
the action ()2.2p is invariant under fohation-preserving diffeomorphism (|2.8p . 
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We point out that the parameter 771 ( 772) is associated with deviation from the BD theory 
along the normal (spatial) direction of the leaves of the foliation. Moreover, we introduced 
two BD parameters ui and U2 which need not to be the same in the anisotropic case. With 
the choice of the parameters of A = 1, ?7i = 772 = 1, and ui = 002 = wbdj the action 
fl2.5p reduces to (isotropic) Brans-Dicke action fl2.2p . For a constant aBD scalar field, it 
becomes Lorentz violating Einstein-Hilbert action with an anisotropic parameter A. Note 
also that the aBD action (12. 5p is still invariant with respect to global (isotropic) conformal 
transformation 

where VL is an arbitrary constant. We further remark that for VL = Q{t,x.), the conformal 
symmetry under the transformation (12. 9p can be extended to local conformal symmetry if 
one chooses 

r/i = -l(l-3A), u;^ = ^{l-3X), 

V2 = 1, C02 = (2.10) 

which implies that for A = 1 case, i.e., 771 = ?72 = 1, oJi = U2 = —3/2, the aBD action (12. 5p 
reduces to a conformally invariant isotropic action [T5] . 

3 Scalar graviton mode in the quadratic action 

In order to investigate the behavior of graviton mode in the Minkowski background {gij = 
6ij, (f) = 0o); we first consider the following scalar perturbations of the metric and aBD 
scalar field for the Minkowski background up to the linear order 

iV = e^ Ni = d,P, gi,=e-^^6,,, <j) = <j)o + V- (3-1) 

In these perturbations, a, /3, ip, and (f are functions of space and time. Substituting the 
above perturbations into the action (12.50 . one finds the quadratic action (□ = did^) 

^(2) = J rftci3^|30o(l - SX)ip^ + 20o(l - SX)ipni3 + 0o(l - A)(n/3)^ + 2r]i{3ipip + ipD^) 
+uji<P^^^'^ + A^Ui) + 40oan^ - 20oV'nV' - 2r/2an^ - a;20o ^(^i^)^}- (3-2) 

Variation of the fields a and (3 of the quadratic action leads to the (local) Hamiltonian and 
momentum constraints 

2 

ip = 0, (3.3) 

up + f 1 - 3A + = 0, (3.4) 

1 - A V r]2 J 
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respectively. 

From the constraints ( 13.31) and (13.41) . we can replace the fields ip and with if). After 
taking the integration by parts, the quadratic action (13. 2p can be rewritten as 

^(2) = 2 j dt£x |-^^' + + ^ ~ ^) ' ^2-^) 

where 

^ 00 I vi V2 vi J 

We notice two things. First is that the perturbation ip can be completely eliminated due 
to the Hamiltonian constraint (13. 3p in favor of the scalar graviton mode. The other is that 
m the hmiti ?72 — > C)0, we find the action is exactly the same with the quadratic action in 
the HL gravity (for 0o = 1) 

SSI = 2 j dtd^xj-^^/'^ - VDV'I , (3.7) 

with c^L = (1 ~ ^)/(3A — !)• It can be easily checked that in the quadratic action (13. 7p . 
the scalar graviton ip has serious problems, i.e., classical instability in the case of c^l < 
or the presence of ghost when c^l > 0. Interestingly, in our framework such pathological 
behaviors do not show up with 

4<0, 0<^ + --l, 3.8 

which can be satisfied with a wide range of the parameters region. 

Let us focus on the special cases rji = 772 (otherwise noticed) which have rather special 
consequences. For rji = r]2 = rj, the quadratic action (13. 5p becomes 

S^''> = 200 I dtd'x { (3 + ^) + + ^ - 1) V^nV^} , (3.9) 

and it shows that the dependency on the parameter A for the scalar graviton completely 
disappears. This also implies the quadratic action ( 13. 5p with (13. 6p is devoid of the singularity 
associated with the limit A — )■ 1. The condition ( 13. 8 p of avoiding the classical instability 
and ghost mode in this case is given by 

coi>-^V^ a;2>-2 + ^(r/-2)2. (3.10) 



^This limit of 7^2 — >■ 00 is reminiscent of a — 00 which is the coefHcient of (WiN/N)'^ term introduced 
in Ref. [16] (see also Ref. [17] for details). 
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Note that uj2 always has to be greater than —2 for any value of rj. For 77 = 1, wi > — | and 
UJ2 > — |, which coincides the ghost free case of the isotropic BD theory [18], r/i = 772 = 
1, ui = U2 = wbd > — f • Another case of interest is when t] = 2. In this case, we have 
Ui > —6, (X'2 > —2. We will see that for this particular value, U2 does not have to be a 
large value in order to pass the solar test as in the original BD theory [see Eq. fl5.7p ]. As r] 
becomes larger, Ui can be more negative, whereas UJ2 shifts away from the value —2 in the 
positive direction. 



4 Strong coupling in the cubic action 

In the previous section, we investigated the behavior of graviton mode in the quadratic 
action and showed that the pathological behavior can be cured in a wide range of parameters 
with the condition of Eq. (13. 8p being satisfied. Another important task is to examine if the 
theory gets strongly coupled at the cubic action level in the limit of A — )■ 1 [T9|[20]. 

In order to study with the cubic order interaction terms [21], we first recall the non- 
linear scalar perturbations (13. ip around the Minkowski background and substitute this into 
the action (12. 5p . After some tedious manipulations one can find that the cubic-order action 
is given by 

^(3) = 

dtd^x{ - 2ip{dipy + 2(f)o'^{d'^y - 2a0o(5^)^ - 4:ipipnip + AaipDip + 20oa^n^ 



-F20oV'^n'^ - 4:(f)oa'ipOip - 'r]2a'^Dip - r]2ip'^D(p + 2r]2ailjDLp - 2r]2ipdiLpdiilj 
+2r]2adi^pdi'ilj + Ws^o V(5<^)^ - (^2<Pta{d^f + ws^o V(5<^)^ - 9(1 - ^X)^^^'^ 
-2(1 - 3A)0o^V'n/3 - 2(1 - ?>\)Md,^dil3 - 2(1 - \)<PoUmW + M{d^djl5f 
+A(j)odidj^di^djij - A0o^(n/3)2 - 3(1 - 3A)0ott^^ - 2(1 - 3A)0oa^n/3 
-a(j)o{didj(3f + (f)oXa{n^f + 3(1 - 3X)ipi'^ + 2(1 - 3A)(pVO/3 + ^{didj^f 



+di^diPDp') - (^(j)o'aip^ + 3(/)o Vvi' + </'o + 2(l)^^ipd,ipdiP'^^. (4.1) 

Using the first-order Hamiltonian and momentum constraints (13. 3p . (13.41) obtained in the 
previous section, the above action (14. ip reduces to 

5(3)=0„/" dtd'x\AM^^^? + A2^^i4^^J^]+A,^(^ + A,^4^'\ . (4.2) 
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where 



^3 

^4 



^2 4^2 12 _ 



r/2 m 



8r/i+4(l-3A)r/2 . 2 . , .ox .^ 2 



(-1 + Xfr^l 
3r]2-2 



2r]f - 4r]ir]2 + (3A - l)r]^ + 2(A - l)wi 



1 



2r/i + (1 - 3A)r/2 



.1 + xy^jli^^ + ^)(3A - l)(3r/2 - 2)7/2^ - 4r/i2(2A + 3r/2(A - 2)^ 

-4r;i772(2-6A + 3r/2(l + A)) - 4wi(A - l)2(3r/2 + 2)}. 



Note that when 772—7-00 the above action (14.21) once again reduces to the cubic action 
in the HL gravity as 



+- 



2 + -^ijd.M ' 



-HL 



□ 



HL 



2chl + 1^^2 



c 



HL 



(4.3) 



It is well-known that the above action is confronted with the strong coupling problem in 
the limit A — >■ 1 [T9||20]. On the other hand, in the case of 771 = 772 = 77 of our main focus, 
the above coefficients A2, A3, /I4 considerably simplify to become 



= 12 3 



2ui 



A, 



-9(3-- 
V 



Aa = 9-- 



6 4a;i(2 + 3r/) 



which again shows that there is no dependency on the parameter A for the scalar graviton 
up to the cubic order. To see this more closely, we introduce the canonical variable ip = 
\/2tjj/\c^ \ with |c^|~"^ = 0o(3+2ci;i/?7^), in terms of which the quadratic action ( 13. 9p becomes 



5(2) 



dt£x <i + ( ^ + ^ - 1^ ^olcv^rv-nv^l 



(4.4) 



By using the canonical variable ijj, the cubic action (14. 2 p can be written as 



5(3) 



2 



dtd^x 



^"^^ + ^ + - - 2) Uc^H{dSf + l2\cS^^^^^ ^ 



-9|3--j0o|c^|'^(^^' 



2 



7/ 



27+- Uo|c^|'-2(3 + - ) |c^| 



□ 



(4.5) 
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This action clearly shows non dependency on the parameter A and that we do not need to 
consider any fine-tuning for the interactions to be regular. 

Before closing this section, we comment on the case with local conformal invariance. It 
turns out that with the choice of the parameters as in fl2.10p . the time dependent cubic 
terms of ip in the action fl4.2p vanish and we get 

5(3) = 00 y dtd^x {-%i){di^l)f] . (4.6) 

This shows that the strong coupling problem does not show up at the cubic-order pertur- 
bation of the action with the local conformal invariance. 



5 Cosmological constraint on anisotropic BD theory 

In this section we investigate the cosmological tests, which can provide the observational 
constraints for alternative theories of gravitation. For this purpose, we first consider the 
solar system test. Taking into account the static part of the perturbations in Sec. 3, the 
Lagrangian (13. 2p with a static point-like source term of mass Mg can be written by 



/^static = {4^nV^ + 4(f)oantlj - 2(f)otljntP - 27]2anip - uj2(j)o\di^y} - M,a5=^(x). (5.1) 
From varying for the fields ijj, a, and (f, we obtain the following equations: 

+ (poDa - (poOijj = 

A(f)onip -2r]2n(p - Ms6^{yi) = 

2n4j - ri2na + uj2(po^a<^ = O. (5.2) 
The corresponding solutions to Eq. (15. 2 p are 

, V2+(^2 (r/2-2)0o .p. ON 
2 + UJ2 002 + 2 

where a is 

a = ("^ + ^)^- (54) 

87r0o(2w2-r/i + 4r/2)|x| 

It is well-known that the linear expansion of the metric (13. ip can be expressed by the 
Newtonian potential u and the Post-Newtonian parameter 7 as follows: 

^00 = -iv2 = _i _ 2a = -1 + 2u, (5.5) 

g,, = il-2ij)6,, = il + 2^u)6,,, (5.6) 
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which imphes that for the solution (15. 3 P u = —a and the parameter 7 can be written by 

7 = — = = 1 H . 5.7 

' a 2 + 002 002 + 2 ^ ^ 

Let us concentrate on the case rji = 1^2 = V again. The solar system tests currently 
constrain 7 as [7] 
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1| < 2 X 10"^ (5.^ 



For the case of the BD theory (77 = 1, 002 = wbd), we see that the above condition (15. 8p 
together with (15.71) restricts the region of oo^d to 

wbd > 50000. (5.9) 

Note, however, that unlike the BD theory, this large value constraint on cubd (15. 9 P can be 
circumvented in the anisotropic BD case when 77 is very close to 2 as can be seen from 

|r/-2| < 2 X 10-^(a;2 + 2) (5.10) 

obtained from Eq. (15.70 . Especially, when 77 = 2, any 002 > —2 is allowed which is also 
consistent with the stability and no-ghost condition of Eq. (I3.10p . In addition, for the 
solution (15. 4p of a one can find the effective Newton's constant Ge// as 

r - ^2 + 2 

being obtained from the Newtonian potential, i.e., u = GeffMs/\x.\. In the General Rela- 
tivity limit {u2 = wbd 00) we have Gefj — G^ for the substitution of 0o ^ l/(167rG7v) 
with the Newton's constant G^. It is interesting to notice that in the anisotropic BD theory 
the reduction Geff G^ can be achieved alternatively with appropriate choices of the free 
parameters different from the BD theory. Especially, when choosing the rj = 2 once again, 
Geff reduces to the Newton's constant G^ irrespective of the value of 002. 

On the other hand, we can also obtain the cosmological constraint from a ratio factor of 
Gc/Geff {Gc being the effective cosmological gravitational constant defined in Eq. (I5.12p ) 
which is related to the primordial helium abundance [23]. Substituting the FRW metric 
ansatz and = (j){t) into the action (12. 5p with matter term included and varying the action 
with respect to the lapse function N[t), we find the standard Friedmann equation can be 
written as 

H' = (5^2) 

where H is the Hubble parameter, p is total matter density of the Universe, and the effective 
cosmological gravitational constant Gc is given by 

Gc = A 7- (5.13) 

87r0(3A-l) ^ ' 
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It is pointed out that Gc is equivalent to Gn for the substitution of = 0o 1/ {^Qt^Gn) 
in the hmit of A — )■ 1. In our case for (p = (pQ, we obtain the cosmological constraint from 
the observational bound [23] of Gc/Geff as 



Gr 



G. 



3(A 



3A 



(3A-l)(u;2 + 2) 



< 0.125. 



(5.14) 



Comparing (15.101) with the bound (I5.14p by using (13.101) . we find that the allowed range of 
the parameter A is 



< A - 1 < 0.095 ^ 10 



-1 



(5.15) 



which imposes a rather loose constraint around A 



6 Conclusion and discussion 

In this paper, we constructed an anisotropic Brans-Dicke gravity, which includes five free 
parameters, i.e.. A, r]i, 772, cui, and uj2. In the case of A = r/i = 772 = 1 and ui = U2 = ojbd, 
the gravity reduces to the ordinary Brans-Dicke action with a BD parameter cjbd- When 
fixing the aBD scalar field to be a constant value, it becomes Lorentz-violating Einstein- 
Hilbert action with an anisotropic parameter A. 

We found that in the perturbation around the Minkowski background and constant aBD 
scalar field the scalar graviton at the quadratic as well as cubic order in the aBD action 
does not show any pathological behaviors within some parameter range. This suggests that 
the aBD gravity can be a viable theoretical model. Especially, the case 771 = 772 reveals 
intriguing property of A independence of both the quadratic and cubic actions. Also in 
the context of cosmological model, we have checked that unlike the case of the BD theory 
which is imposed by cubd > 50000 to be consistent with the experimental observations, this 
large value can be evaded in the aBD theory with the special value of 7/ = 2 which is one of 
the novel feature of the aBD theory. But the origin of these special properties associated 
with 7/ = 2 is a puzzling aspect which needs to be investigated further. 

We conclude with a couple of comments on the issues related to the BD theory. First, 
it is well-known that in the standard BD model without potential, there is no accelerated 
expansion, so one has to consider a potential term [9l[10]. However, it is found in the aBD 
model including matter contribution that we have a de Sitter solution {H = const.) for the 
FRW metric, given by 

a = aoe^\ = 0oe^=^(i+--)^*, = poa^'^'^^-^ (6.1) 

with constants oq, 00, Po and Pm = oUmPm- In particular, A and po satisfy the following 
relation: 

X = 27]+{1- cuDcui + 1/3, Po = -9um(l)oaf^^"'^H\7] + uj,{1 + u^)). (6.2) 
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For Pm = (vacuum) case, the above relation yields 

= -2r/ + A - 1/3. (6.3) 

It should be pointed out that in the BD limit, i.e., ?7 — ?■ 1, A — )■ 1 the parameter ui = ubd 
becomes a negative value, ubd = —4/3 which conflicts with the lower bound ubd > 50000, 
even if it satisfies the ghost-free condition of ubd > —3/2 [18]. However, in the aBD case 
this problem can be circumvented by choosing ?7tobe?7>4/3or?7<0 for the bound 
(I5.15p . obtained from substituting (16. 3p into the ghost-free condition of Eq. (I3.10p of Sec. 
3, when rji = 772. It is worth noticing that the allowed range is not in conflict with the 
special case of 77 = 2. 

The second one is a speculation about quantum gravity. In Refs. [25] - [27] . the one-loop 
counter terms in the pure BD theory is calculated and it is shown that the BD theory is 
not renormalizable unless curvature squared terms are included. Likewise, we suspect that 
adding curvature squared terms [l2l[T6l[2Hl[29] in the aBD theory might constitute the UV 
completion of the theory. The details are beyond the scope of the present paper, but the 
subject deserves further investigations. 
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